Adjustable subwavelength localization in a hybrid plasmonic waveguide 
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The hybrid plasmonic waveguide consists of a high-permittivity dielectric nanofiber embedded in a low- 
permittivity dielectric near a metal surface. This architecture is considered as one of the most perspective 
candidates for long-range subwavelength guiding. We present qualitative analysis and numerical results which 
reveal advantages of the special waveguide design when dielectric constant of the cylinder is greater than the 
absolute value of the dielectric constant of the metal. In this case the arbitrary subwavelength mode size can 
be achieved by controlling the gap width. Our qualitative analysis is based on consideration of sandwich-like 
conductor-gap-dielectric system. The numerical solution is obtained by expansion of the hybrid plasmonic 
mode over single cylinder modes and the surface plasmon-polariton modes of the metal screen and matching 
the boundary conditions. 

PACS numbers: 



I. INTRODUCTION 

The creation of the waveguides capable of guiding light 
with deep subwavelength confinement is of great interest for 
practical applications. These devices may throw open the 
doors to nanoscale optical communications, quantum comput- 
ing, nanoscale lasers and bio-medical sensing. The main prob- 
lem on the way to practical realization is the diffraction limit 
of light in dielectric media. Electromagnetic energy cannot be 
localized into nanoscale region much smaller than the wave- 
length of light in the dielectric 1 1 1. The possible solution to 
this problem is using of the materials with negative dielectric 
permittivity. For example, metals are known to exhibit this 
property below the plasma frequency. Metal structures pro- 
vide guiding of the surface plasmon-polaritons (SPP), which 
can be strongly localized near metal-dielectric interfaces (3. 
However the propagation length of the strongly confined plas- 
monic modes is not large enough due to the presence of Ohmic 
losses in the dissipative metal regions. 

The new approach for this challenge integrates dielec- 
tric waveguide with plasmonic one. The hybrid plasmonic 
waveguide consists of a high-permittivity dielectric nanofiber 
separated from a metal screen by low-permittivity dielec- 
tric nanoscale gap |3|. Both the single fiber and the silver- 
dielectric interface cannot provide strong mode confine- 
ment at optical and near infrared frequencies, but such hy- 
brid conductor-gap-dielectric architecture has experimentally 
demonstrated deep subwavelength optical waveguiding |4|. 
Relatively large propagation distance has been achieved due to 
low loss tangent tg = s'^/s'^ at the operating frequency and 
the specific spatial structure of the guiding mode with field 
confinement within non dissipative gap region. 

In the present paper we show that the hybrid plasmon po- 
lariton(HPP) mode confinement can be considerably risen by 
a specific choice of the materials, when the dielectric constant 
of the cylinder is greater than the absolute value of the di- 
electric constant of the metal screen. The main advantage of 
the choice is the hyperbolic-like dependence of the effective 
index on the gap width. This feature allows to achieve ar- 
bitrary subwavelength mode size at any frequency by tuning 
the distance between the cylinder and the metal. To justify 



our approach we theoretically investigate propagation of the 
HPP-mode. First we give qualitative analysis basing on the 
consideration of plane sandwich-like conductor-gap-dielectric 
waveguide structure(CGD) |5 1. We derive exact analytical ex- 
pression for effective index of the fundamental CGD-mode 
and give criterion when the HPP-mode is CGD-like. Fi- 
nally, we present semi-analytical approach for describing of 
the HPP-mode propagation. Similar approach has been previ- 
ously applied for plane wave scattering by a cylinder placed 
near the plane surface tSl-fTll. 

The scheme is based on expansion of the HPP-mode over 
single cylinder modes and the surface plasmon polariton 
modes of the metal-dielectric interface and matching the 
boundary conditions for electromagnetic field components. 
Numerically obtained dispersion relations confirm the advan- 
tages of our design of the hybrid waveguide. 
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FIG. 1: a) Geometry of the waveguide; b) Plain waveguide with 
the same width of the gap. 



II. QUALITATIVE DESCRIPTION 

The geometry of the hybrid waveguide is the following: a 
circular dielectric cylinder of diameter d and permittivity £d 
is placed above a metal screen of permittivity Sm- The width 
of the gap between the cylinder and the metal screen is h, 
see Fig. [TJa). Let us choose the Cartesian reference system 
as it is shown in Fig. [T} z-axis is directed along the waveg- 
uide, whereas x-axis is directed normally to the metal screen. 
We consider a plasmon-polariton mode of frequency uj and 
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the propagation constant /3 propagating along z-axis. Thus, 
all electromagnetic field components depend on time and z- 
coordinate as exp[z/3z — iuji\. We assume, that responses 
of both dielectric and metal on electromagnetic field are de- 
scribed by dielectric constants, which are and respec- 
tively. Generally, the outer medium may be not vacuum, but 
some dielectric medium having dielectric constant being equal 
to Sg. All the materials are assumed to be nonmagnetic. To 
describe the mode confinement, it is convenient to introduce 
effective refractive index neff, which is defined as n^^ — ^jk, 
where k = uj/c'm the wavenumber in vacuum. The effective 
index determines the field penetration depth into the material 
with permittivity e as 1 //c \/^eff ~ ^- penetration depth of 
the bound mode should be real in the unbounded waveguide 
constituents (metal and outer dielectric space), and may be 
imaginary for bounded constituents (fiber). The greater neff is 
the stronger degree of confinement. 

Optimization for transversal field confinement imple- 
mented in paper 1 3 1 for hybrid waveguide shows that the thin- 
ner gaps provide higher localization. The fiber diameter is 
much greater than the the gap width in the case, and the mode 
is sufficiently localized in the region where the gap can be con- 
sidered as approximately plain. In the region, the waveguide 
shape is close to plain sandwich like conductor-gap-dielectric 
(CGD) structure, see Fig.[Tj6). The limit of plain CGD-model 
was considered in 0, where the properties of the bound fun- 
damental mode were investigated. One of the main advan- 
tage of the CGD-mode is that effective index of the mode 
ncGD is greater than the refractive index of the dielectric 
^CGD > ^^^^ implies, that the electromagnetic field 

of the mode decays exponentially into the dielectric cladding. 
Nevertheless, the analysis proposed in |5 1 is not applicable to 
HPP-mode of hybrid waveguide with optimal diameter found 
in p|. The reason is that the mode of plain CGD-model in- 
deed describes the HPP-mode only for large enough fiber di- 
ameter d otherwise HPP-mode should be considered as a re- 
sult of hybridization of surface plasmon polariton modes and 
the modes of the single dielectric cylinder. 
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The main goal of the present work is to give the theoretical 
description of the hybrid waveguide and to find approaches 
to deeper localization of the HPP-mode. Comparative analy- 
sis of O] [5l suggests, that in order to get stronger transversal 
miniaturization of the hybrid waveguide the CGD-like regime 
of propagation (with neff > -s/sd) should be achieved for the 
diameter which is much less than free space wavelength. Our 
analysis of the plain CGD-structure shows that the localiza- 
tion of the fundamental mode can be significantly risen for 
special set of materials, when absolute value of the metal di- 
electric constant is less than the dielectric constant of the di- 
electric cladding, < £d- For the case, the effective refrac- 
tive index ncoD is proportional to inverse width of the gap, 
ncGD 1/kh, when the width h is small enough. To use 
the same effect for the hybrid waveguide, the cylinder diam- 
eter should sufficiently exceed some critical value d*, which 
is determined by the condition that the transversal size of the 
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FIG. 2: (Color online) Effective index ticgd of the CGD-mode ver- 
sus gap width h for different metal permittivity and critical gap width 
{syn] he). The dielectric constants of the dielectric and gap region are 
Sd = 5.76 and £g = 1 respectively at wavelength A = 490nm. 



plain part of the gap is comparabl e with th e mode penetra- 
tion depth into the dielectric l/k^/nlf^ — Ed- The size of the 



plain part of the gap is evalua ted as 2 v hd, thus the condition 
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■Ed- For d greater enough than 



the guiding mode can approach the strongly confined mode of 
the sandwich like system even if the diameter of the cylinder 
is much less than free space wavelength. 

In order to give general physical argumentation of our re- 
sults let us consider planar sandwich-like CGD-waveguide in 
detail. The wave vector of fundamental CGD-mode (which is 
TM-mode) can be calculated from equation | 5 1 



ex.p[2hf<ig 



{SdK.g + Sgl^d){Sml^g + SgKrn) ' 



(2) 



where hii = k^nl^^^ — Si for each material, i = m^g^d and 
^cGD is the effective index of the mode. In particular, 1/hZd 
and I /him are the penetration depths into the dielectric and the 
metal correspondingly. It is known that such plane three-layer 
waveguide supports the propagation of the bound eigen mode 
only if the width of the intermediate layer is less than some 
cut-off value he which is determined by the permittivities at 
given frequency 

When the thickness exceeds this critical value the funda- 
mental CGD-mode becomes radiative with energy leaking 
into upper dielectric half space. 

There is a significant difference between the dependence of 
the effective index nccD on the gap thickness h for the cases 
of low and high index dielectric, see Fig|2] For relatively low 
refractive index of dielectric, Ed < \£m\, there exists surface 
plasmon-polariton mode when the gap is absent, h = 0. It 
has effective index nmd = \/^m^d/ (^d + ^m)- Then the 
effective index of the fundamental CGD-mode is bounded, 
^/£d < ^cGD < ^md- Just the case was considered in the 
papers |5| and |3|. Otherwise when permittivity of dielectric 
is relatively high 



^d > km > ^g 



(3) 
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the effective index ncoD unHmitedly diverges as the gap thick- 
ness tends to zero, h <C A/y^: 



The diameter dependence of the effective index 



1 (£d -gg)(gm -gg) 

2kh {ed + gg)(gm + £g) ' 



(4) 



This leads to extremely strong light confinement in a trans- 
parent dielectric gap layer located between the high-index di- 
electric and the conductor. The actual degree of localization is 
restricted only by additional factors, such as increasing Ohmic 
losses in the metal, spatial dispersion and atomic structure of 
the materials. In this respect the properties of the conductor- 
gap-dielectric plasmonic mode similar to that of the gap plas- 
mon polaritons in a conductor-gap-conductor structure |8|- 
191 . This feature is the principle behind our idea: in prac- 
tice one should choose the metal of the absolute permittivity 
less than the permittivity of the cylinder and place cylinder 
at distance h < he from the metal plane. When such metal 
is involved the effective index of the HPP mode can be sig- 
nificantly greater than effective index of electromagnetic field 
in bulk material of the cylinder even for very small diameters 
of the cylinder. Note, that to calculate group velocity Vg and 
chromatic dispersion for the mode using the formula (|4]), one 
should know the dispersion laws for permittivities Em and ^d- 
For thin gap h <^ the group velocity scales 8iSVg/c(xh/X. 
Thus divergence of the CGD-mode effective index with gap 
width decreasing leads to strong reduction of the group veloc- 
ity. 

There is reverse side of the strong localization which is 
small propagation distance. It was shown in paper |3| that 
the strongest localization of the HPP-mode corresponds to the 
lowest propagation length. It is common place of waveguides 
which use metal as a constructive component. Let us consider 
limit when the gap index is low, so Eg <C |£m | • For the case 
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where e^^ is the imaginary part of the metal permittivity. It 
follows from Eq. ([5]), that the localization radius is of the or- 
der of h\em\ in the limit h ^ he. Note, that our approach 
allows to squeeze the mode at arbitrary frequency into any 
subwavelength scale simply by tuning the gap width in accor- 
dance with (|4]). Hence, our waveguide design breaks connec- 
tion between mode localization and the carrying frequency of 
the mode. In particular, the approach may be interesting for 
design waveguides at THz frequencies ITOl |TT1. The prop- 
agation length £ ~ h\em\/\tg\ i.e. reduces with the mode 
size reduction. To keep the propagation length acceptable for 
practical implementation at fixed degree of localization one 
should minimize loss tangent tg. Thus, a prospecting like \TT\ 
is needed to propose the optimal choice of materials for our 
approach ([3]). 



III. SEMI-ANALYTICAL DESCRIPTION 

In the section, we present the semi-analytical approach to 
the propagation of the HPP-mode and discuss the numerical 
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FIG. 3: (Color online) Effective refractive index of the fundamental 
hybrid mode versus cylinder diameter d (coloured lines) compared 
with those of single fiber (black solid line) and SPP mode (lower black 
broken line). The dielectric constants of the cylinder, dielectric and 
metal are Sd = 12.25 Sg = 2.25 and Sm — —129 + 3.3i respec- 
tively at wavelength A = 1.55/im. These parameters are chosen in 
accordance with the paper |3 |. The critical gap width he — 5nm. 
The HPP-to-CGD crossover point: ^ 17/im for h = 2nm. 



results. It follows from Maxwell's equations that the elec- 
tromagnetic field of guiding mode can be fully described by 
z-components of the electric and the magnetic fields, and 
Bz |13|. These fields satisfy the following two-dimensional 
Helmholtz differential equation inside the homogeneous areas 
where permittivity is constant: 



= 0, 



(6) 



where A-^ 



dl + and k 



uj/c is the free space 

wavenumber. The boundary conditions on the both interfaces 
are continuity of components Ez, Hz, sE^ and H^, where ^- 
component of a vector is its normal component. 

Our semi-analytical method is based on the representation 
of the hybrid waveguide as an integration of the dielectric 
fiber and plane plasmonic waveguide. We express the elec- 
tromagnetic field of the HPP-mode as a linear combination 
of cylindrical modes around the fiber and evanescent plane 
waves above the metal screen. Boundary conditions provide 
the system of linear equations on the expansion coefficients. 
Such an approach leads to highly efficient method of numer- 
ical solving a difficult boundary- value problems that describe 
the propagation of waves in a complex systems |[T4ll - |[T5l . The 
scheme is developed in detail in Appendix [A| 

To verify our semi-analytical method, in Fig. [3] we present 
the dependence of the effective index of the fundamental hy- 
brid mode on the cylinder diameter d for a range of the gap 
widths h in the case of telecommunication wavelength when 
< < km|. These dispersion curves are obtained from 
our numerical procedure and show a good agreement with the 
results obtained in |3J by using finite-element package FEM- 
Lab from COMSOL. 
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The diameter dependence of the effective index 
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FIG. 4: Wave vector of the fundamental hybrid mode versus cy Under 
diameter d (coloured lines) compared with those of single fiber (black 
solid line) and SPP mode(lower black broken line). The dielectric 
constants of the cylinder, dielectric and metal are £d = 5.76, Sg = 1 
and Sm = —9.2 respectively at wavelength A = 0.49/im. These 
parameters are chosen in accordance with the paper |16 1. The critical 
gap width he = 7nm. The HPP-to-CGD crossover points are: ^ 
SlOnm for h — 2nm , ^ 875nm for h — 5nm. 



FIG. 5: Wave vector of the fundamental hybrid mode versus cylinder 
diameter ^(coloured lines) compared with those of single fiber (black 
solid line) and SPP mode(lower black broken line). The dielectric 
constants of the cylinder, dielectric and metal are = 5.76, Sg = 1 
and Sm = —4 respectively at wavelength A = 0.49/im. The critical 
gap width he = 13, 4nm. The HPP-to-CGD crossover points(black 
arrows) are: ^ 40nm for h — 2nm, (f ^ 65nm for h — 5nm, 
(T ^ 220nm for h = lOnm. 



In accordance with general argumentation given in Sec- 
tion[ll|we next present two sets of plots. Fig.Qa) corresponds 
to the case of fiber with comparatively low refractive index, 

< km|, the parameters of the waveguide are taken accord- 
ingly to experimental work 1 16|. Fig.|4j6) corresponds to op- 
posite limit, when > km|- Parameters of these two plots 
differ only for metal permittivity e^, the value = —4 is 
chosen for Fig.|4]^6). Here, we do not concretize the material 
of the metal screen, our goal is just to demonstrate the qual- 
itative difference of the guiding mode properties for the case 

Presented results indicate that when fiber diameter d is de- 
creased, the HPP-mode loses confinement along the metal and 
eventually (at (i = 0) becomes a surface plasmon-polariton 
mode of the flat metal- vacuum interface. Herewith the ef- 
fective index of the HPP-mode monotonically decreases to 
that of this SPP-mode. Thus all dispersion curves have the 
same asymptotic neff n^g = \/^m%/ (^m + at small 
d. Two different behavior are possible at the opposite limit 
of large diameter. As the diameter d oo, the HPP-mode 
can asymptotically tend either fundamental single fiber mode 
or the fundamental mode of the planar three-layer system, 
the choice depends on the gap width h. If the gap thick- 
ness h is below than h^^ (Eq.([T])) the HPP-mode approaches 
the CGD-mode with the diameter increasing. In the case the 
crossover between the asymptotics occurs at (i* (black arrows 
at Fig.lJF^)) which is determined as 



6^* 



1 



(7) 



For h > he the HPP-mode becomes the cylinder-like 
in the limit of the large diameter. In the case the critical 
diameter do corresponding to the transition between small- 



diameter and large-diameter asymptotics is defined by the 
equation nsp((io) = ^mg, where ngp((i) is the diameter de- 
pendence of the effective index of the single fiber fundamen- 
tal mode. If the condition ^/e^kd <C 1 is valid one can de- 
rive that the localization of this mode is exponentially small. 



exp 



eg{ed-e^){kdY 



<1, (8) 



and 7 = 0.5772... is Euler-Mascheroni constant. 

Let us suppose the effective index nccD of the CGD-mode 
not to be significantly above that of bulk plane wave in the 
fiber medium, y/e^. For example, this assumption is true at 
the conventional plasmonic condition when the absolute value 
of metal permittivity is sufficiently greater than the dielec- 
tric permittivity, Jsm I /^d — 1^1- Just the case is realised 
at Figjsjand Fig|4ja). Then the field penetration depth into 
the upper dielectric is quite large as well as the HPP-to-CGD 
crossover diameter d* > A, so the CGD-mode does not pro- 
vide the strong confinement. Therefore there are no advan- 
tages of CGD-like limit in the case from the view of HPP- 
mode confinement. For a given frequency and gap width the 
choice with the strongest coupling of the fiber mode and the 
surface plasmon polariton mode, corresponding to d = do, 
provides the strongest localization of the field within nanogap 
due to the great contrast of permittivities |3 |. At the same 
time significant part of energy is transferred inside the fiber, 
thus the waveguide mode confinement is achieved largely due 
to the boundedness of the high-permittivity dielectric part of 
the waveguide. Once the diameter of the fiber is optimum 
and the gap width is small enough the advantages of the hy- 
brid architecture are used completely: cross section size of 
the system can be much less than the wavelength and mode 
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confinement is much stronger than for uncoupled single fiber 
or flat metal-dielectric interface. To achieve further increase 
of the HPP-mode confinement the fiber with higher dielectric 
constant should be used. 

Next let us assume that effective index of the CGD-mode 
is significantly larger than the refractive index of the fiber 
medium. This can be achieved by diminishing the gap thick- 
ness in the case \£m\ < Sd that corresponds to the Fig|4]^6). 
Then the CGD-mode has strong confinement so the crossover 
diameter can be decreased to deep subwavelength scale, d* <C 
A, by tuning the gap width. Therefore the attractive CGD-like 
asymptotic can be achieved by HPP-mode with very small 
diameter of cylinder providing the wished structure of the 
mode with the strong transversal localization in two dimen- 
sions within the gap region and exponential decaying into the 
cylinder. Note that in the case the top part of the fiber cross 
section is at distances much larger than from the gap and 
its particular shape does not play role any more. 



IV. CONCLUSION 

In the paper we have proposed the novel approach for hy- 
brid plasmonic waveguide design providing wide opportuni- 
ties for HPP-mode property controlling. When the absolute 
permittivity of the metal is less than that of the dielectric the 
hybrid effective index is unlimitedly diverges (Eq.Q) with 
gap width decreasing. High effective index provides strong 
confinement of the electromagnetic field in two dimensions 
within the nanometer- scale gap region. Thus the mode size 
can be simply controlled by tuning the waveguide' geometry 
at fixed frequency and materials constituting the waveguide. 
The advantages of the case \em\ < £d are confirmed by both 
qualitative analysis within planar three-layer model and rigor- 
ous semi-analytical method describing the HPP-mode propa- 
gation in general. The propagation distance of hybrid mode 
reduces with the mode size reduction. To achieve long-range 
propagation at fixed degree of localization one should mini- 
mize loss tangent tg = {e'^/^ml the metal. It should be 
noted that simultaneous satisfying of both conditions \em\ < 
Ed and tg <C 1 at optical and near infrared frequencies is a 
challenging task. Thus implementation of the waveguide loss 
compensation techniques would be required to use such hy- 
brid waveguide as a component of the miniaturized photonic 
circuits. Another potential application of our waveguide de- 
sign lies in study field of the resonant plasmons. The reso- 
nance condition for a surface plasmon-polariton at a planar 
metal-dielectric interface is the fine-tuning of the permittivi- 
ties, — {Sm + Sd) <C Ed iflTl . Thus for particular metal the 
resonance can be achieved only in a narrow spectral range. 
While resonant increasing of the CGD-mode effective index 
requires only the geometrical condition <C A. Thus for 
the conductor-gap-dielectric structure the resonance of plas- 
monic mode can be attained by gap width decreasing at any 
frequency as long as the condition \Em\ < Sd 'i^ valid. 




FIG. 6: Reference system. 



Appendix A: Numerical method 

The theoretical description of the hybrid waveguide is in- 
hibited by its complex geometry. In general we should chose 
such system of coordinates where the surfaces of the waveg- 
uide are the isolines and Helmholtz equation can be solved by 
separation of variables. The hybrid geometry corresponds to 
the so called bipolar coordinates based on two sets of orthog- 
onal circles. In this coordinate system the Helmholtz equation 
has quite complicated form and accordingly the set of eigen 
functions cannot be found analytically. However the unknown 
hybrid eigen functions can be expressed in terms of known 
solutions of the Helmholtz equation in other coordinate sys- 
tems. It is convenient to represent the total electromagnetic 
field of HPP-mode as the superposition of the all modes of sin- 
gle fiber(cylindrical functions) and all SPP modes(evanescent 
plane waves) with some unknown coefficients of expansion. 

We chose the Cartesian system of coordinates with origin 
at axis of the cylinder which is z-axis (Fig. |A]). Supposing 
the structure of the fundamental hybrid mode to be symmetric 
with respect to x-axis we can describe the longitudinal com- 
ponent of the electric field as 

, Ei'^^= EZo b^Kniv) COS n(^+ 

+ /o^ exp(Q/^g(x - D)) cosqn^y dq 
^ Ei"^^ = ex.p{-QK.rn{x - D)) cosqK.rny dq 

(Al) 

and z-components of the magnetic field as 

H^z"^^ = Er=o^n^n(Xdr)sinn(/p 

+ /o cx.p{Qtvg{x — D)) sinqtvgif dq 

H^^^ df ex.p{-QK.^{x - D)) smqK.^y dq 

(A2) 

where Q \J\^q^ and D = d/2^ h. 

To write the corresponding equations, it is convenient to ex- 
press the field inside dielectric in terms of only plane evanes- 
cent waves, when we impose the continuity conditions on the 
boundary of the metal, and in terms of angular harmonics, 
for the cylindrical surface. We solve it by using the evanes- 
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cent plane wave expansion of modified cyHndrical functions 
and angular harmonic spectrum of the evanescent plane waves 

m 



Kn{n^r) cos flip 



/ F^{q)e-'^-^^cosq^^ydq, (A3) 
Jo 



Kn{f^gr)smn(p = / F^{q)e~^'^^'^ sin qi^^y dq, (A4) 
Jo 



^Qi^gx ^^^y ^ ^ 

n=0 

oo 

^QKgX ^^^y ^ ^ 



(A5) 
(A6) 



n=0 



where 



.E^ {Q + qr + {Q-qr 

2Q 

,H _ (Q + g)" - (Q - g)" 

2Q 



(A7) 
(A8) 



2 - (5o,i 



HQ + qr + {Q - qDUv), (A9) 



((Q + 'Zr-(Q-9)")/nKr'). 



(AlO) 



Thus the electromagnetic fields in surrounding medium 
close to the dielectric waveguide can be written as 

oo 

E^^^ = ^6^i^n(/^gr)cosn(/p + 

n=0 ■'0 



n=0 



n=0 



POO 

+ ^sinn(/;/ c^G^{r,q)e-^^^''dq. 

n=0 

The corresponding expressions for the fields close to the sur- 
face of the metal are 

= / E b^Fniq)e-'^''''' cos 9«gy + 

■^0 n=0 

/>CX) 

+ / c^e^''^^''-^^cosqKgydq, 
Jo 



poo 

^i'^ = / E ^n^n (^)e-^^^" sin^^g^ dq + 

n=0 

' s'mqtZgy dq. 



It can be easily derived from the Maxwell's equations that 
for the normal components 



ik dH, 



13'^ — ek?' d(, — ek'^ dr] 

ik dE^ ip dH^ 
'^(3^-€k^~d^ ~ (3^ - ek^ ~df 



(All) 



, (A12) 



where rj is tangent to the interface coordinate in the transversal 
plane. 

The continuity conditions on the metal surface for E^, Bz, 
eEx and eB^ lead to the first system of linear homogeneous 
equations (SLE) on coefficients 6^, b^, c^, c^, d^, d^. The 
corresponding continuity conditions for Ez, Bz, sE^ and sB^ 
on the cylindrical surface produce the second SLE on ampli- 
tudes a^, a^, b^, b^, c^, c^, which is now integral with re- 
spect to c^, c^. In order to avoid integration of the unknown 
functions we express the coefficients c^, in the terms of 6^, 

from the first SLE and substitute them into the second SLE. 
The procedure leads to the infinite system of linear homoge- 
neous algebraic equations for coefficients a^, a^, b^, b^. In 
order to solve the system numerically one should truncate it 
to a finite size. Then the propagation constant of the funda- 
mental hybrid mode can be determined from the condition of 
vanishing of the characteristic determinant. 
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